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ABSTRACT 

We calculate a grid of spherically symmetric OB stellar atmospheres at low metal- 
licities, including both non-local thermodynamic equilibrium (NLTE) and metal line 
blanketing effects. This is done to assess the uncertainties in helium abundance deter- 
mination by nebular codes due to input stellar atmosphere models. The more sophis- 
ticated stellar atmosphere models we use can differ from LTE models by as much as 
40 percent in the ratio of He to H-ionizing photons. 
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1 INTRODUCTION 

Giant extra-galactic H II regions (GEHR) are interesting 
objects for many reasons (cf. Shields 1990 for a review). 
However, low metallicity HII regions draw the most atten- 
tion. This subgroup of low metallicity GEHR is believed 
to have had little chemical evolution over a Hubble time. 
Therefore, it is widely accepted that GEHR offer an oppor- 
tunity to measure the primordial chemical composition of 
the universe. 

The primordial composition of the universe is predicted 
by big-bang nucleosynthesis theory. If low metallicity GEHR 
provide a way to determine the primordial element abun- 
dance then they allow testing of some big-bang cosmology 
predictions, in particular the primordial helium abundance 
Y p , which has implications for the baryonic fraction of mat- 
ter in the universe and may limit the number of exotic light 
particle species. Although low metallicity GEHR show little 
chemical evolution, some chemical evolution did take place 
over a Hubble time, thereby requiring an extrapolation of 
the "measured" helium abundance Y obs (as a function of 
the observed metallicity) to zero metallicity. This procedure 
yields the value of Y p , the quantity of interest. The extrap- 
olation techniques to zero metallicity are discussed by Pegal 
et. al. (1992). In this paper we concentrate on systematic er- 
rors in the determination of Y ob ° from low metallicity GEHR 
introduced by systematic changes of the radiation field that 
ionizes the nebula. 

The use of nebular spectroscopy to measure element 



abundance is model dependent. The internal physics of 
GEHR is therefore of prime importance. The uncertainties 
include: a) the spectrum of the ionizing low metallicity stars, 
b) the filamentary structure of GEHR, c) supersonic turbu- 
lence inferred from radio observations, d) radiative transfer 
effects within the nebula (Terlevich, Skillman & Terlevich 
1995 ;Sasselov & Goldwirth 1995). 

The uncertainty of Y obs determinations is required to 
be no more than several percent in order to derive accurate 
enough values of Y p for, say, establishing the number of 
light neutrino species. Typically the value of Y° bs is derived 
f rom recombination lin es of singly and doubly ionized helium 
( |3killman et al. 1994 ) which are compared to H a emission. 
This implies that an essential assumption is that the fraction 
of neutral helium in the nebula is negligible. If the fraction 
of neutral Helium is several percent, then the value of Y° bs 
cannot be determined to the required accuracy. Therefore, 
the helium ionization degree in H II regions has been the 
subject of numerous meticulous studies. Two problems are 
often mentioned in this context: i) the filamentary structure 
of the nebula and ii) the uncertainties in the underlying 
ionizing (stellar) radiation field. The latter uncertainties can 
translate into uncertainties in relative sizes of the He and H 
ionization zones. 

This paper does not attempt to deal with problems as- 
sociated with the filamentary structure of the nebula. How- 
ever, the existence of this problem must be noted. It is ob- 
servationally established that most of the line emission of 
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GEHR comes from a rather small filamentary volume. The 
volume filling factor of these filaments is about 1%. They 
move relative to each other with a velocity dispersion which 
is supersonic. Moreover, radio observations reveal that in 
about half of the H II regions in M51 (not a low metallicity 
region) the radio emission originates in a non-thermal com- 
ponent within the filaments. The origin and dynamics of the 
filaments are not fully understood. Clearly, if they contain 
any neutral or singly ionized He, then the He/H abundance 
determination is hindered by this source of inaccuracy. Thus, 
before any final statement can be made, one requires strong 
observational constraints of the properties of the matter in 
these filaments or a reliable theoretical model to understand 
the origin and the detailed physics of the filaments. 

Due to the complexity of the nebular filamentary struc- 
ture most authors have invoked homogeneous nebular mod- 
els, in slab or spherical geometry. These spherical and shell 
homogeneous models constrain the ionization degree of he- 
lium to be 0.98 if the e ffective temperature of the ionizing 



stars is above 38,000K (Dinerstein & Shields 1986 



Skillman 



198S()~]l owever, these authors note that if the GEHR is ac- 
tually an ensemble of much smaller H II regions which are 
ionized by stars of different effective temperatures (lower 
than 38,000K, Dinerstein 1989), filamentary neutral helium 
structures might form within the nebula. Even if this prob- 
lem is ignored and one assumes that all ionizing stars have ef- 
fective temperatures above 38,000K, there are uncertainties 
due to a lack of a state-of-the-art g rid of model a tmospheres 
in the required range of interest (Shields 1990). Recently. 



some e ffort h as been put into this proble m by (Garacia- 
Vargas 1996; Garacia- Vargas et. al. 1997). They attempt 



to account for the filamentation by using shell models of 
t he nebula, th en, using the photo-ionization code CLOUDY 
(Ferland 1996), they fit the st ellar and nebular spectra self - 
consistently. They conclude ( Garacia- Vargas et. al. 1997) 
that contamination from young evolved stellar population 
e.g. Wolf-Rayet stars, is of importance and must be taken 
into account. 

The need for co nsistent spectral modeling of GEHR has 
been recognized by ( Garacia- Vargas & Diaz 1994), who have 



used (Kurucz 1992) LTE line blanketed models with a range 



of abundances as input to CLOUDY. However, they express 
their concern that NLTE effects might become important 
in the relevant range of parameters. This lack of state-of- 
the-art NLTE line blanketed model stellar atmospheres is 
troubling, because most of the ionizing spectra of the nebula 
is in the UV range, which is strongly metal line blanketed. 
Typically for solar abundance stars, photons in this spectral 
region are strongly affected by the presence of many spectral 
lines. Yet, the input stellar model atmospheres used so far 
assumed either a pure hydrogen helium mixture in NLTE, or 
assumed LTE for the metal lines. This introduced a source of 



unc erta inty that this paper addresses. In addition, ( Kurucz 
f 994 ^Jn is improved and modified the line data set of ( Kurucz 
1992), and these effects are taken into account as well. 

There is already a substantial literature on non-LTE 
effects (some of the classical papers are: Auer & Miha- 
las 1972, Kudritzki 1973, 1976, 1979, 1988, Husfeld et al. 
1984). Our purpose here is to evaluate their contribution 
to the uncertainties in He abundance determination. In re- 
cent years stellar atmosphere and radiative transfer codes 
have become very sophisticated in a way which allows much 



more of the important physics to be taken into account. 
We use such a code to construct a grid of low-metallicity 
hot stellar atmosphere models including NLTE and metal 
line blanketing. This allows us to assess and reduce one of 
the uncertainties in the determination of Y° bs . We provide 
the spectra of these models in machine readable format at 



http : //dilbert .physast .uga. edu/~yeti. 

The structure of this paper is as follows: in § ^| we discuss 
briefly the stellar atmosphere code PHOENIX , and describe in 
detail the input physics of stellar spectrum models. We then 
provide in § |^ some examples of the variation of the spec- 
trum in the UV range. We present the ionization parameters 
of helium and hydrogen, depending on various assumptions, 
and assess the variation of these quantities as functions of 
the stellar atmosphere model input parameters. We summa- 
rize our conclusions in § ^. 



2 METHODS AND MODELS 

In order to investigate the importance of NLTE effects on 
the formation of OB star spectra and ionizing photon fluxes, 
full and very detailed NLTE model calculations are required 
in order to model the effects of NLTE on the very important 
EUV and UV metal lines. This means that the multi-level 
NLTE rate equations must be solved self-consistently and 
simultaneously with the radiative transfer and energy equa- 
tions, including the effects of line blanketing and the ex- 
tension of the atmosphere. For the purpose of this analysis 
we use our multi-purp ose stellar atmosphere code PHOEMIX 
PHOENIX [ve r sion 9.1 (|Hauschildt et 



et al. 1997 



al. 1996; flauschildt 



Hauschildt, Baron, 



Allard 1997T Baron & 



Hauschildt 1998 )7 

uses a special relativistic spherical radia- 



tive transfer and an equation of state (EOS) including more 
than 300 atoms and ions (39 elements with up to 26 ioniza- 
tion stages). The temperature correction is based on a mod- 
ified (for NLTE and scattering) Unsold-Lucy m ethod that 
converges quic kly and is numerically stable. See ( Hauschildt 
& Baron 1996 ) for details on the numerical methods used in 
PHOENIX. 

Both the NLTE and LTE (background) lines are treated 
with a direct opacity sampling method. However, we do not 
use pre-computed opacity sampling tables. We dynamically 
select the relevant LTE background lines from master line 
lists at the beginning of each iteration and sum up the con- 
tribution of every line to compute the total line opacity at 
arbitrary wavelength points. The latter is crucial in NLTE 
calculations in which the wavelength grid is both irregu- 
lar and variable (from iteration to iteration due to changes 
in the physical conditions). This approach also allows de- 
tailed and depth dependent line profiles to be used during 
the iterations. To make this method computationally effi- 
cient, we employ modern numerical techniques, e.g., vector- 
i zed and parallel block algorithms with high data locality 
(Hauschildt, Baron, & Allard 1997), and we use high-end 



workstations or parallel supercomputers for the model cal- 
culations. In the calculations we present in this paper, we 
have set the statistical (treated as micro-turbulence) veloc- 
ity £ to 2kms~ . We include LTE background lines (i.e., 
lines of species that are not treated explicitly in NLTE) if 
they are stronger than a threshold F = Xi/ K c = 10 -4 , where 
Xi is the extinction coefficient of the line at the line center 
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and k c is the local b-f absorption coefficient. This typically 
leads to about 2 x 10 6 LTE background lines. The line pro- 
files of these lines are taken to be depth-dependent profiles 
(with Voigt profiles for the strong lines and Gauss profiles 
for weak lines). We have verified in test calculations that the 
details of the LTE background line profiles and the thresh- 
old r do not have a significant effect on the model structure 
and the synthetic spectra. However, the LTE background 
lines should be included in detailed model calculations be- 
cause their cumulative effect can change the structure and 
the synthetic spectra. In addition, we include about 2000 
photo-ionization cross sections for atoms and ions (Verner 
& Yakovlev 1995). 

PHOENIX is a full multi-level NLTE code, i.e., NLTE ef- 
fects are considered self-consistently in the model calcula- 
tions, including the radiative transfer calculations and the 
temperature corrections. Hauschildt & Baron (1995) have 
extended the numerical method developed by Hauschildt 
(1993) for NLTE calculations with a very detailed model 
atom of Fe II. In the calculations presented in this paper, 
we use a significantly enlarged set of NLTE species, namely 
H, He I-II, Mg II, Ca II, Ne I, C HV, N I— V I, O I- VI, S II- 
III, and Si II — III, (Hauschildt et al. 1997, for a complete 
list of NLTE species available in PHOENIX 9.1 see). Here, 
we do not use the NLTE treatment for Li I, Na I, Fe I— 
III, Co I — III and Ti I — III because these ionization stages 
are not important in metal-poor OB star atmospheres and 
particularly Fe, Co, and Ti NLTE would considerably in- 
crease the CPU time for the model calculations with little 
additional improvements in the results. We include a total 
of 2120 NLTE levels and 14,080 NLTE primary lines in the 
calculations presented here, nearly a factor of 5 more levels 



scribed in (|Hauschildt 1992"|). Details of the method are de- 
scribed in (Hauschildt & Baron 1995), so we do not re- 



and line s than in our previous calculations (Hauschildt et al 
1996). The construction of the model atoms is described in 



(Hauschildt et al. 1997) and references therein. 



2.1 NLTE Calculation Method 

The large number of transitions that have to be included 
in realistic models of the NLTE line formation require an 
efficient method for the numerical solution of the multi- 
level NLTE radiative transfer and model calculation prob- 
lem. The model atoms used here include more than 14,080 
individual NLTE lines plus a large number of weak back- 
ground transitions. Classical techniques, such as the com- 
plete linearization or the Equivalent Two Level Atom meth- 
ods, are computationally prohibitive. PHOENIX performs its 
calculations in the co-moving frame for expanding atmo- 
spheres (e.g., stellar winds, novae and supernovae), there- 
f ore, approaches such as Ande rson's multi-group scheme 
( Anderson 1987 ; Anderson 1989 ) or extensions of the opacity 
distribution function method (Hubeny & Lanz 1995) can- 



not be applied. Simple approximations such as the Sobolev 
method, are very inaccurate in problems in which lines over- 
lap strongly and make a significant continuum contribution 
(important for weak lines), as is the case for nova (and SN ) 



atmospheres cf. (Hauschildt et al. 1996; Baron et al. 1996) 
We use the multi-level rate-op erator splitting (o r pre- 
conditioning) method described by ( Hauschildt 1993 ). This 
method solves the non-grey, spherically symmetric, special 
relativistic equation of radiative transfer in the co-moving 
(Lagrangian) frame using the operator splitting method de- 



peat the detailed description here. For all primary NLTE 
l ines the radiative rates and the approximate rate operators 
(Hauschildt 1993) are computed and included in the itera- 
tion process. Secondary NLTE lines are included as back- 
ground transitions for completeness but are insignificant for 
the model structure and the synthetic spectra (the model 
atoms have been explicitly constructed so that all impor- 
tant lines are primary lines). 



2.2 Atmosphere Models 

We have computed a small grid of spherically symmetric 
static models to investigate the effects of NLTE on the struc- 
ture and the spectra of metal poor OB-stars. The models in- 
clude the NLTE treatment as discussed above as well as the 
standard PHOENIX NLTE generalized equation of state and 
additional LTE background lines (about 2 million atomic 
lines). Aufdenberg et. al. (1998) showed that the combined 
effects of line blanketing and spherical geometry can dramat- 
ically change the short wavelength spectrum of the mod- 
els, therefore, we calculated all models presented here us- 
ing spherical geometry (incl. spherically symmetric radia- 
tive transfer). Models that include the effects of a stellar 
wind on the structure of the atmosphere and the emitted 
spectrum are currently being calculated (Aufdenberg et al, 
in preparation). The NLTE effects are included in both the 
temperature iterations (so that the structure of the mod- 
els includes NLTE effects) and all radiative transfer calcu- 
lations. For each primary NLTE line we add 3 to 5 wave- 
length points to the overall wavelength grid. This procedure 
typically leads to about 55,000-150,000 wavelength points 
for the model computations and the synthetic spectrum cal- 
culations. Test calculations have shown that the resulting 
overall wavelength grid is completely adequ ate for NLTE 
calculations in static and expanding media ( [Hauschildt & 
Baron 199^ ). We have also calculated LTE continuum and 
line blanketed models for comparison. In figure |l| we show 
a synthetic spectrum at the nominal resolution, illustrating 
the effect of low resolution on the appearance of the spec- 
trum. 

We compare our LTE and NLTE models to the Ku- 
rucz 92 LTE set of synthetic spectra in Figs. |^ and ^. The 
PHOENIX spectra have been convolved with a Gaussian ker- 
nel of 6A half-width to make the resolutions of the different 
sets of spectra comparable. In general the spectra are very 
similar. We use an updated version of the Kurucz atomic 
line lists; this could explain the differences between the LTE 
models. In addition, the PHOENIX models allow for line scat- 
tering and NLTE (Fig. |) whereas the Kurucz models are 
calculated using complete LTE (no line scattering). 



3 RESULTS 

All model atmospheres presented in this paper have log(<?) = 
4.0. Three effective temperatures (A, B, and C respectively) 
T ef f = 38, 000, 45,000, and 55, 000 K, have been considered. 
For each of these T e g we have computed models with 10%, 
5% and 2% solar metal abundances(l-3, 4-6, 7-9 respec- 
tively). In each case NLTE and LTE line blanketed mod- 
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els as well as LTE continuum models were calculated. The 
models and the resulting ionization parameters go: the base 
ten-logarithm of the emission rate of HI Lyman continuum 
photons 



10 = logio i 4tt" 



dX 



he 



(1) 



and qi (helium continuum) are presented in table 1. In addi- 
tion we have computed for each temperature a line blanketed 
LTE solar metallicity model (AS, BS, CS), to assess quan- 
titative differences which could result from a contaminating 
population of of young evolved stars that may be present in 
the HII region. 

We start by considering low metallicity models. As ex- 
pected, the resulting ionization photon fluxes are most sensi- 
tive to the temperature. The dependence in gi on the metal- 
licity is typically weak. The strongest variation is among the 
different model type themselves (LTE vs. NLTE). For the 
higher temperatures the differences in go and gi are typi- 
cally of order 0.02 or less. For models Al-9, the differences 
in gi between LTE and NLTE models can be as large as 
0.19, meaning that the amount of He ionizing flux differs by 
a factor of 10 019 . As the differences in q is only about 0.02, 
the implied difference is the ratio of helium-ionizing flux to 



hydrogen-ionizing flux is about 10 



1.45. 



We compare the spectra of the NLTE, LTE line blan- 
keted and LTE continuum models for three effective tem- 
peratures in Figs. ^ to ^. The effects of line blanketing on 
the overall shape of the spectra is significant, even at the 
lowest metallicity that we have considered. The differences 
between LTE and NLTE spectra are smaller than the effects 
of line blanketing by itself. Longward of 1000 A, NLTE ef- 
fects on the spectra are small, especially near the threshold 
temperature of T e « ~ 38, 000 and would thus have small if 
any impact on near-UV observations. 

In figures Q to |^ we show comparisons between solar 
metallicity spectra and spectra computed for 1/10 of the 
solar metallicity. In the solar metallicity spectra the metal 
lines are always stronger than in the spectra for 1/10 solar 
metallicity (though the Lyman edge is typically not quite 
as strong as in the 1/10 solar metallicity spectra). The dif- 
ferences in the value of gi are 11%. This implies that if 
10% of the stars are younger higher metallicity stars, the 
differences in the total ionizing number of photons excit- 
ing the surrounding nebula would be about 1%. However, 
the differences between the spectra are considerably more 
pronounced. The re sulting effects from this fac t have been 
studied recently by ( Saracia- Vargas et. al. 1997), who found 



that the Wolf-Rayet contribution typically dominates this 
effect. 



4 CONCLUSIONS 

Our synthetic spectra show that the predicted degree of he- 
lium ionization in GEHR could vary significantly depend- 
ing on different model assumptions about the stellar atmo- 
spheres. In the end, the computed value of [He]/[H] in the 
nebula depends on the fraction of ionization contributed by 
the "relatively cool" hot stars, as well as the geometry of 
the HII region. As the number of recombinations in an HII 
region equals the number of ionization, one might expect 



that the depths of the ionization zones for He and H would 
scale roughly as the cube root of the ionizing fluxes (only 
roughly because He-ionizing radiation can also ionize H) . On 
the other hand, the uncertainties discussed in this paper, i.e. 
those due to the assumption of a static atmosphere, might 
be augmented by additional uncertainties in line blanketing 
effects, such as the existence of winds (Gabler et al., 1989, 
1991, 1992; Najarro et al. 1996), shocks, etc., which are be- 
yond the scope of this paper. The existence of additional un- 
certainties in the nebular dynamics, such as those implied 
by the presence of filaments, also need to be investigated 
farther. The extent to which the theoretical uncertainty can 
be gauged by scatter in the data is also hard to assess here, 
because much of it may be systematic. 
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Table 1. Models Input Parameters and Results. Cont stands for 
continuum LTE models, LTE for line blanketed LTE model and 
NLTE for NLTE line blanketed models 



Model 


Type 


T eff 


[M/H] 


90 


91 


Al 


NLTE 


38,000K 


-1. 


24.0443 


23.1560 


A2 


LTE 


38,000K 


-1. 


24.0232 


22.9831 


A3 


Cont 


38.000K 


-1. 


24.0226 


23.0756 


A4 


NLTE 


38,000K 


-1.3 


24.0351 


23.1610 


A5 


LTE 


38.000K 


-1.3 


24.0204 


22.9928 


A6 


Cont 


38.000K 


-1.3 


24.0307 


23.0976 


A7 


NLTE 


38,000K 


-1.7 


24.0311 


23.1885 


A8 


LTE 


38,000K 


-1.7 


24.0164 


22.9959 


A9 


Cont 


38.000K 


-1.7 


24.0292 


23.0991 


AS 


LTE 


38.000K 





24.0640 


22.9930 


Bl 


NLTE 


45,000K 


-1. 


24.5376 


23.8949 


B2 


LTE 


45.000K 


-1. 


24.5345 


23.9063 


B3 


Cont 


45,000K 


-1. 


24.5113 


24.0245 


B4 


NLTE 


45.000K 


-1.3 


24.5322 


23.9171 


B5 


LTE 


45.000K 


-1.3 


24.5345 


23.9063 


B6 


Cont 


45.000K 


-1.3 


24.5100 


24.0246 


B7 


NLTE 


45.000K 


-1.7 


24.5282 


23.9497 


B8 


LTE 


45.000K 


-1.7 


24.5247 


23.9512 


B9 


Cont 


45.000K 


-1.7 


24.5092 


24.0246 


BS 


LTE 


45,000K 





24.5620 


23.818 
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Figure 1. The NLTE spectrum for T eS = 45, 000 K, log(g) = 4.0 
and solar metallicity at the nominal resolution used in the model 
calculations. 
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Figure 2. Comparison between a PHOENIX NLTE spectrum and a 
Kurucz 1992 LTE spectrum for T cff = 45, 000 K, log(g) = 5.0 and 
1/10 solar metallicity. The PHOENIX sspectrum has been degraded 
in resolution by convolving it with a Gaussian kernel with 6A 
half-width. 



Figure 3. Comparison between a PHOENIX NLTE spectrum and 
a Kurucz 1992 LTE spectrum for T off = 45, 000 K, log(g) = 5.0 
and solar metallicity. The PHOENIX sspectrum has been degraded 
in resolution by convolving it with a Gaussian kernel with 6A 
half-width. 
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Figure 4. Models with T eff = 38kK, Continuum LTE, Line 
LTE and NLTE models are overlaid. Top panel 10% solar metal 
abundance, bottom panel 2% solar metal abundance. 
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Figure 5. Models with T e ff = 45kK, Continuum LTE, Line 
LTE and NLTE models are overlaid. Top panel 10% solar metal 
abundance, bottom panel 2% solar metal abundance. 
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Figure 6. Models with T e ff = 55kK, Continuum LTE, Line Figure 9. Comparison between Synthetic spectra for T e ff = 

LTE and NLTE models arc overlaid. Top panel 10% solar metal 55kK for solar and 2% solar metal abundances, 

abundance, bottom panel 2% solar metal abundance. 
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Figure 7. Comparison between Synthetic spectra for T e f j = 
38kK for solar and 2% solar metal abundances. 
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Figure 8. Comparison between Synthetic spectra for T e ff = 
45kK for solar and 2% solar metal abundances. 
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